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Abstract 

We study effects of imposing various forms of the kinematical constraints on the full form 
of the CCFM equation and its non-linear extension. We find, that imposing the constraint 
in its complete form modifies significantly the shape of gluon density as compared to forms 
of the constraint used in numerical calculations and phenomenological applications. In 
particular the resulting gluon density is suppressed for large values of hard scale related 
parameter and kx of gluon. This result might be important in description of jet correlations 
at Large Hadron Collider within CCFM approach. 


Introduction 

In the high energy limit of hadron scattering, in the regime where the center of mass energy is 
larger than any other available scale, perturbative approach to processes with high momentum 
transfer allows factorization of the cross section into a hard matrix element with initial off-shell 
gluons and an unintegrated gluon density m- The unintegrated gluon density is a function of 
the longitudinal momentum fraction x and transverse momentum kx of a gluon. After taking 
into account formally subleading corrections coming from coherence of gluon emissions, one 
is lead to the CCFM set of equations EH5] which introduce gluon density dependent on hard 
scale related to the probe as well as one introduces unintegrated quark densities. The status of 
phenomenological relevance of CCFM framework is not fully satisfactory. In principle, it is a set 
of equations that should be the ideal framework for application to final states at high energies 
and covering DGLAP and BFKL domains. It has been implemented in the Monte Carlo event 
generator [6]. However, so far good agreement with high precision data has been successfully 
achieved only in rather inclusive processes like F-i and Drell-Yan [7]. It is known that on the 
theory side the CCFM physics is still to be completed. Below we outline the main points. For 
instance: 

• the CCFM system of equations has been so far solved in decoupled approximation in 
the, gluonic (for short the equation for gluon density is usually called CCFM) and non¬ 
singlet channels neglecting correlations between partons [8|. This might result in improper 
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Figure 1: The plot visualizing kinematical variables used in the CCFM equation. The black blob 
represents hard process while the gray oval represents hadron. 


treatment of the kinematical region where the gluon degrees of freedom are less dominant 
than sea quarks. Since gluon might be artificially dominating over quarks where it should 
actually be suppressed. Furthermore the corrections coming from inclusion of transversal 
momentum dependence to splitting functions P gq and P qq are not known and in order to 
have complete picture they should be calculated [§]. 

• the impact of kinematical effects introducing energy conservation in the CCFM evolution 
have not been investigated in all detail. As it turns out from our study it is necessary to 
revisit the inclusion of the so called kinematical constraint [9] into the CCFM equation. 
The orgin of the kinematical constraint follows from refining the assumption that the t- 
channel off-shell gluon’s 4-momentum is dominated by k^- The kinematical constraint 
causes suppression of the gluon density and even overrides the angular ordering in regions 
of large kx- 

• only recently the CCFM has been promoted to non-linear equation |10ffT2] therefore al¬ 
lowing for the possibility to investigate interplay of coherence effects and saturation [1| in 
exclusive processes like dijet production at the LHC [13U20| . In particular the important 
question is what is the role of the angular ordering and kinematical effect in the evolution 
at the non-linear level. The optimal form of initial conditions is also not known. 

This publication is a continuation of the work done in |21jj22]. We compare numerically forms 
of non-Sudakov form factor used in the literature and solve the full CCFM equation and its 
non-linear extension - including the 1/(1 — z) pole and kinematical constraint in the kernel of 
the equation i.e. not only in the non-Sudakov form factor. However we keep a s constant in 
order to have clear picture of the role of the kinematical effects. 


The CCFM equation 


The CCFM equation for gluon density reads: 


A(x, k T ,p) = Ao(x, k T ,p ) + a s 


d 2 q 
7rq 2 



9(p - zq ) V(z, k T , q ) A s (p, zq , Q 0 ) A(x/z, k T q), 

(1) 


where k' = |k + (1 — z)q| and the moduli of the two dimensional vectors transversal to the 
collision plane are denoted |k| = kx, |q| = qx, x is gluon’s longitudinal momentum fraction and 
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a s = N c a s /ir. Also the rescaled momentum is introduced as q = |q| = qy/(l — z). 
The splitting function assumes the form: 

-r,/ , — \ A NS (z,k^,q) 1 

V{z,k T ,q) = ---+ 


2 1 — z 

The Sudakov form-factor evaluated in double logarithmic approximation reads: 


( 2 ) 


A s (p,zq,Q 0 ) = exp - 


dq 2 

q ' 2 


i Q o 
|q'l 


dz! ■ 


«5 

1 -z’ 


(3) 


' ( Z q)2 0 

The non-Sudakov form-factor regularizing 1 jz singularity in angular ordered region is: 


Aivs(z, hr,q) = exp ( “ ^ 


1 , kn 


= exp ( —as In - In 


T 

z zq 1 


(4) 


The kinematical constraint 

The integration over q in the equation 0, although being constrained from below by the soft cut¬ 
off Qo, is not constrained by an upper limit thus violating the energy-momentum conservation. 
Moreover in the low x formalism one requires that in the denominator of the off-shell gluon 
propagator one keeps terms that obey \k 2 \ = fey. In order to be consistent as argued in [3j the 
non-Sudakov form-factor should be accompanied by a kinematical constraint limiting the above 
integration over q. In approximated form it reads 


fey > zq 2 


(5) 


After including it in the (J4]) we have: 


Ans(z, k T , q) = exp 



( 6 ) 


The condition Q at z <S 1 guaranties that \k 2 \ ~ k^. In [9j it has been extended to region 
including also the case when z ~ 1. Below we re-obtain the full form of kinematical constraint 
emphasizing the its role in conservation of energy. Having k = z p + + zp~ + k± with p + and p~ 
being the initial state gluon momenta + and — components, we can write the expression for k 2 


k 2 = —z z s — kj' , 


(7) 


with s = ( p + + p~) 2 = 2 p + ■ p~. Note, that for the full propagator \k 2 \ > fey, thus the full 
propagator causes stronger suppression of the amplitude which has to be taken into account 
when using \k 2 \ ~ fey. A cut-off like (0 is a simple implementation of the suppression. 

The condition \k 2 \ ~ fey translates approximately to 

fey > z z s. (8) 


Using the identity q 2 = z( 1 — z)s — qx 2 = 0, we can express zs = qT 2 /{ 1 — z ) and insert it 
into © to obtain: 

fey > j = z (1 — z) q 2 (9) 

The lower bound on z > x results in the upper bound on qtf, < fcy/x ~ s providing local condition 
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Figure 2: The plots visualizes the various forms of non-Sudakov form factor. On the left we choose the 
situation when kx = 10 GeV and q = 4 GeV while on the right kx = 4 GeV and q = 10 GeV. The red 
line corresponds to form-factor given by formula the brown line corresponds to form-factor given 
by formula (ED while blue line to form-factor given by formula (1121) . The form-factor given by ED is 
multiplied by the function 6 (kx/ (1 — z)q 2 — z) 


for energy-momentum conservation. The CCFM with the kinematical constraint included reads 

M- 


A(x, k Tl p ) = A 0 (x, k T ,p) + a s 


[cfti f 

J Trq 2 J 


i Q o 
1 |q| 


dz 6 


k^ 


— z 


.(! - z)q 2 

X 

x 0(p-z\q\)V(z,kr,q) A. s (p,zq,Q 0 )A{x/z,k T ',q) , 


( 10 ) 


and the non-Sudakov form-factor after inclusion of the full form of the kinematical constraint 
assumes form 


&ns(z, k T ,q) = exp - a s 


f 1 dz ' 
f dq’ 2 

J q ’ 2 


0 


(1 — z')kx 


(1 — z) 2 q 2 


— z 


Q{k 2 T - q' 2 ) Q(q' - z'q) 


( 11 ) 


Please note the presence of the function 9 ^ ~ z ^j • At the level of the kernel of the 

BFKL equation obtained directly from Feynmann diagrams (virtual and real corrections are 
kept apart) one imposes the kinematical constraint at the so called "unresummed" level on real 
emissions. However in |9j it has been observed, at a level of BFKL equation, that imposing 
the kinematical constraint on unresummed equation and performing algebraic transformations 
leading to resummed equation one obtains a 0-function affecting the full kernel. The starting 
point of the CCFM is a resummed kernel, therefore this result suggested how the complete form 
of the kinematical constraint should be imposed on the kernel of CCFM. The authors of [9] solve 
the CCFM at small z limit therefore the function 6 — z ^j is neglected and in most of the 

phenomenological and theoretical applications of the CCFM this term is neglected [7l ll5i[23H28j 
as well. The following form of non-Sudakov form factor is usually being used: 

A Ns(z,kT,q ) = exp (j-as J ^ J ^®(k 2 - q ,2 )Q(q'- z'q)^j (12) 

= exp f -a s log ) log 



4 
















where 


f 1 , if (kr/q) > 1 

zo= l k T /q, if 2 < (kr/q) < 1 (13) 

[ z , if {kr/q) < z 

and the discussed above the ^-function is not taken into account. 

As shown in the Fig. [2] the simplified non-Sudakov form-factor © approximates well the exact 
one (EH) in regions of z where the 0-function is not affecting the kernel while the (1121) is quite 
different. In our analysis we consider following equations: 


A(x, k T ,p) = Ao(x,k T ,p) + a s J J dz 0 ^ _^- 2 “ z 

X 


X 

x 9(p — zq) V(z, k T , q) A s (p, zq, Q 0 )A(x/z,k T ',q) 
where the splitting function includes non-Sudakov given by formula (ED and 


A(x,k T ,p) = A 0 (x,k T ,p) + a s [ 

J vrq 2 


i_l2a 

Q 


7rq z 

X 

V d(n — yn\ / D( v bm m 


dz 


(14) 


(15) 


where the non-Sudakov given by formula (1121) . 


Saturation effects and kinematical constraint combined 

To account for gluon recombination at large gluon densities the CCFM equation has been pro¬ 
moted to non-linear equation by including a quadratic term |10Wl2| which reads: 


i_£a 
<? 


A(x,kr,p) = A 0 {x,k T ,p) + a s J j dz 6 - z\6(p - zq) V(z,k T ,q) 


X 

/ -\ 


(16) 


X &s(p,zq,Q o) [ A(x/z,k T ',q) - 5 \q~- 


k T 2 


(i- z y 


q 2 A 2 (x/z,q,q ) 


where we also included the kinematical constraint of the form (ED in the kernel. Simpler versions 
of the equation above have been already analyzed in |22| and it has been observed that 


the equation leads to phenomenon called saturation at the saturation scale [221126] 
the saturation strongly suppresses the gluon density at low x and low hr 


The natural question arises how are these results modified when some of the approximations are 
not taken and how are they modified if the kinematical effect is imposed in the full form. 


Numerical results and the discussion 

We present the numerical results with as = 0.2 and the parameter Qq set to Q o = 1 GeV. To 
study the variation of the results depending on Q o we produce a solution also for Q o = 0.5 GeV. 
The choices of the starting scale are motivated by a possibility to compare the equations in a 
region where both linear and non-linear equations are phenomenologically relevant. 
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Figure 3: Red (dashed) - equation (fl4l) . black (solid) - equation (11511 . blue (dotted) - equa¬ 
tion (flfil) . 


p=50 GeV, x=1CT 4 /c,=10 GeV, x=1CT 4 




Figure 4: Red (dashed) - equation (fT41) . black (solid) - equation (1151) . blue (dotted) - equa¬ 
tion (flfil) . 


Linear equations 

We use an initial condition which includes resummed virtual and unresolved contributions, 
according to [22j and m , in the form 

Ao(x,kr,p) = AA R (x,kx,Qo) A s (p,Qo,Qo)/kx, ( 17 ) 

where A = 1/2 and A R {z, kx,Q o) = exp ( — as log 4 log -^z J is the Regge form-factor. 

The first observation we make from the plots in Fig. [3] [5] is that the solutions of equations we 
study differ significantly. The solutions exhibit also similar features. Solutions of both versions 
of the kernel with kinematical constraint exhibit a local maximum as functions of kx and p 
with x and p or kx correspondingly fixed. The positions of local maximum in the plots of p 
dependence are correlated with the value of kx, with a shift to higher kx for the solution of 
the equation (11511 . The peak can be explained by the fact that the contribution of the integral 
on the right hand side of CHD peaks at around kx ~ p ■ To point out: the peak is a result of 
presence of 9 (p — zq ) - angular ordering condition. Similar peaks are present also in the plots 
of k R dependence and resemble Sudakov suppression of kx scales of the order of p in case of 
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Figure 5: Red (dashed) - equation (fl4l) . black 
tion (flfil) . 


V=3 GeV, p=3 GeV 



Figure 6: Red (dashed) - equation (fl4l) . black 
tion (flfil) . 



(solid) - equation (11511 . blue (dotted) - equa- 


V= 10 GeV, p=3 GeV 



(solid) - equation (USD , blue (dotted) - equa- 


equation (H3D I203. However, in the case of (USD it seems, that the position of the peak does 
not depend on the value of p. It seems that the kernel with the kinematical constraint included 
in an implicit way (1151) produces similar effect, though weaker, as an implicit inclusion of the 
0-function only in the case of p dependence, but for the for dependence leads to much smaller 
suppression. Thus the peak observed in solution of (USD is ’buried’ under the result of the 
evolution. We can conclude that the peak in the fop dependence is a result of an interplay of the 
inclusion of the explicit 6 ^ 2 — factor and the Sudakov effect. 

Non-linear equation 

We set the parameter characterizing the strength of the non-linear term R the value R = 
\J 1 /7r GeV in the equation (1161) . 

By comparing the CCFM and KGBJS equations we see that the equations give quite similar 
distributions. This effect (for our choice of the starting scale Q 0 ) comes from the fact that 
the kinematical constraint suppresses the growth of the gluon so much that the non-linear 
effects enter only at very low x. Observations made in previous paragraphs are confirmed in 
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Figure 7: Relative ratio of CCFM (1141) and KGBJS (1161) solutions. Distributions with definite 
kx for varying value of p. 


p=10 GeV p=20 GeV 



kj_ [GeV] k i [GeV] 

Figure 8: Relative ratio of CCFM (fT4j) and KGBJS (fl6l) solutions. Distributions with definite p 
for varying value of kx- 


2-dimensional plots (Fig. [THE]), where we plot absolute relative difference of two amplitudes, 
solutions of the CCFM and the KGBJS equations, defined by the quantity 


P(x,k T ,p) 


\Aqcfm (x, kx,p ) — Akgbjs (x, kx,p) 
■Accfm (^, kx,p ) 


(18) 


The function (3 ( x , kx,p ), introduced before in [22], can be used to measure the strength of the 
non-linear effects and to define a saturation scale using the condition: 


(3 ( x,Q s ( x,p ) : p) = const. 


(19) 


or p-related saturation scale: 

/3 ( x , kx, P s ) = const. (20) 

The conditions above can be seen as equipotential lines in the 2-dimensional plots in Fig. [7] [8] 
where different equipotential lines correspond to different constants on the right-hand side of 
the equation above. The change in the slope of the /3 (x, Q s (x,p) ,p) at around kx = p reported 
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p=10 GeV, x=1CT 4 p=50 GeV, x=10' 4 



Figure 9: Solutions compared with Qo = 0.5 GeV and Qo = 1 GeV. Black solid - linear CCFM 
with Qo = 0.5 GeV, red dashed - non-linear CCFM with Qq = 0.5 GeV, blue dotted - linear 
CCFM with Qq = 1 GeV, brown dashed longer - linear CCFM with Qo = 1 GeV. 



P [GeV] p [QeV] 


Figure 10 : Solutions compared with Qo = 0-5 GeV and Qo = 1 GeV. Black solid - linear CCFM 
with Qo = 0.5 GeV, red dashed - non-linear CCFM with Qo = 0.5 GeV, blue dotted - linear 
CCFM with Qo = 1 GeV, brown dashed longer - linear CCFM with Qo = 1 GeV. 


in [22] , apparent in Fig. Ena can be understood m the context of the peuk ut p r ^ j k'j 1 (Fi^. |3F 
ED. For kx > p the contribution of the integral on the right-hand side of the CCFM equation 
decreases and the gluon density is dominated by the initial condition. 

By comparing the plots Fig. [T] [8] to analogous plots in } 22j we see that their main features 
are very similar. We therefore conclude that the low — x approximation of the KGBJS and 
CCFM equations taken in [[22] does not, at least, modify the relative difference between linear 
and non-linear equation since the modifications of the kernel did not spoil the saturation pattern 
visible in the full equation. 

Dependence on the starting scale 

The starting scale Qo dependence is studied in Fig. l9lfT0l using equations f[T4l) and (1161) . We 
observe variation from few up to 25% mostly in the local maximum of the kx and p distributions. 
There is a peculiarity of the dependence on Qo of the kx distributions near the boundary kx = Qo 
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which can be explained by the dependence on the initial condition on the parameter Q q. Near 
the point kx = Q o the initial condition dominates and since it depends on the parameter the 
solutions also depend strongly on this parameter. We can tell, though, that for kx 3> p and 
kx the difference diminishes. 

Conclusions and outlook 

We have solved linear and non-linear version of the CCFM equation with the full kernel, constant 
a s and kinematical constraint included in two different ways. 

Our numerical results make us conclude that the inclusion of the kinematical constraint in its full 
form into the CCFM equation seems to affect largely the distribution of gluons. In particular 
inclusion of the kinematical constraint without omitting the 0-function in the kernel of the 
equation causes strong suppression of the solution for all values of transversal momenta kx- The 
effect is stronger than the implicit kinematical constraint included via modification of the non- 
Sudakov form-factor widely used in literature. We expect that the suppression of the distribution 
will have significant effects on description of decorelations in azimuthal angle of forward-central 
jets [151129] . In particular we expect that the inclusion of the 0-function is crucial for suppression 
of the cross section of emission of gluons in small angles since this configuration corresponds to 
large values of kx- However, the conclusive comparison of the different approaches presented in 
this article can be done only after fitting the initial conditions of the considered equations to 
data, this we, however, leave for future tasks. 
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